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Introducing a perturbative definition, phase space path integrals can be calculated 
without slicing. This leads to a short-time expansion of the quantum-mechanical 
path amplitude, or a high-temperature expansion of the unnormalized density ma- 
trix, respectively. We use the proposed formalism to calculate the effective classical 
Hamiltonian for the harmonic oscillator. 



1 Introduction 



^ I Path integrals are usually evaluated by time-slicing , since the continuum 

definition is mathematically problematic. This becomes obvious for physical 
systems with nontrivial metric, where reparametrization invariance seems to 
be violated. 

An alternative method of calculation was proposed on the basis of a per- 
turbative definition of path integrals in configuration spaced motivated by 
^^ , perturbative procedures in quantum field theory. In this approach, the path 

C^ ' integral of any system is expanded about the exactly known solution for the 

free particle in powers of the coupling constant of the potential. 
C*^ ' The extension of this definition to functional integrals in phase space is 

^^ , presented in this paper by treating the complete Hamilton function as per- 

turbation, which leads to a short-time expansion of the quantum-mechanical 
path amplitude. In Euclidean space, the density matrix is obtained as a high- 
temperature expansion. By simple resummation, this series can be turned 
►v> , into an expansion in powers of the coupling constant of the potential de- 

^ ■ scribed above. As will be shown in the sequel, the knowledge of an exactly 

5^ , known nontrivial path integral as that of a free particle is, however, not re- 

quired. Thus, the perturbative definition presented here is more general. It 
reproduces the expansion around the free particle by a simple re-summation. 
The method is then applied to calculate the effective classical Hamiltonian 
of the harmonic oscillator Hi^^csiPQ,xo) by exactly summing up the pertur- 
bation series. This quantity is related to the quantum statistical partition 
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function via the classical looking phase space integral 

dxodpo 



2-Kh 



exp {-/3iJ„,cff(po, a;o)} , 



(1) 



where /? — l/ksT is the inverse thermal energy. 



2 Perturbative Definition of the Path Integral for the Density 
Matrix 

Slicing the interval [0,h(3] into A^ + 1 pieces of width e = h(3/{N + 1), the 
unnormalized density matrix can be expressed in the continuum limit by^ 



N 


r /"^^ 1 


N+1 


11 


/ dXn 


11 


n=l 


\_J — oo ] 


n=l 



dpn 

2T:h 



Jp„(Xn-Xn-l) 



: exp S -e ^ H{pn, Xn)/h \ , 



(2) 



where Xa — xo and Xb — xn+i are the fixed end points of the path. Expanding 
the last exponential in powers of s/h, we recognize that the zeroth order 
contribution to the density matrix ^ is an infinite product of S -functions due 
to the identity 



dpn 

2TTh' 



,ip„(a;„-ri;„_i) _ 



5{xn -a;„_i). 



(3) 



This infinite product simply reduces to 



lim 



dXN ■ ■ ■ dxi d{xN+l — Xat) • • • 6{X2 — Xl)S{xi —Xo)= S{xb — Xa)- (4) 



Note that this is the unperturbed contribution to the unnormalized density 
matrix ([2]) which was obtained without solving a nontrivial path integral. 
Thus, the phase space path integral for the unnormalized density matrix ^ 
can be perturbatively defined as 



g{xb, Xa) = S{xb - Xa) + y2 ^n ■ / dvi ■ ■ ■ dr^ 

X (H(p(Ti),a;(Ti))..-i/(p(r„),x(r„)))o (5) 
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with expectation values 

N r 

(•••)o 



lim TT / 



dXr, 



N+l r 

n 

n=l 



dpn 
2TTh 



Jp„{x„-x„-i) 



(6) 



These expectation values are usually reexpressed by Feynman diagrams. This 
is possible for polynomial as well as nonpolynomial functions of momentum 

o 

and positioEf^. 



3 Restricted Partition Function and Two-Point Correlations 

The trace over the unnormalized density matrix ([S]) of our unperturbed sys- 
tem with vanishing Hamiltonian H(j), x) = leads to a diverging partition 
function. It is also obvious that the classical partition function diverges with 
the phase space volume. The regularization of these divergences is possible by 
excluding from the phase space path integral the zero frequency fluctuations 
xo and po of the Fourier deco mposition of the periodic path x{t) and mo- 
mentum p(r), respectively^'J^. Therefore we express the quantum statistical 
partition function of any system by 



Z 



dxodpa 
2Trh 



ZPo 



(7) 



where the restricted partition function was introduced as the Boltzmann factor 
of the effective classical Hamiltonian: 

ZP«"« = exp{-/3iJeff(po,a^o)} = 27r?i i VxVpS{xo -x)6{po -p) 



X exp ■ 



hl3 



dr 



'i{p{T) 



-Po)^(a;(r) 



Xa) + H{p(T),x{T)) 



with the measure 



N+l r 



jVxVp^Yu.^\{ J 



dXndpn 



2TTh 



,(8) 



(9) 



nhfl 



and temporal mean values a; = L dTx{T)/'h(3 a.ndp = L dTp{T)/'hl3. As in 
the preceding section, the unperturbed system may have a vanishing Hamil- 
tonian, i.e. H = 0. The calculation of the restricted partition function Zq°^° 
for this system, i.e. evaluating the path integral ^ ior H = 0, is in the same 
sense trivial as for its density matrix, since it also reduces to a cancellation of 
(5- functions and yields Zq"^" = 1. 
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In what follows, we concentrate on the correlation functions of position 
and momentum dependent quantitites. For this purpose it is convenient to 
introduce the generating functional 

ZP°''° [j, v] = 2nn <j> VxVpS{xo ~ x)5{po - p) 



X exp ■ 



np 



dr 



-ip{T)—x{T) +j{t)x{t) +v{t)p{t) 



(10) 



with abbreviations x(t) — x{t) — xq and p{t) = p{t) — pq. The calculation 
yields 



Zrij,«]=exp|^^ 



dr dT'j{T)GP'"'°{T,T')v{T') 



(11) 



where the periodic Green function is 



GP°'-'°{t,t') 



2/3 



{2(t - r') - hl3 [e(T ~ t') - e(T' - r)]} 



2i ■^-^ sina;m(T — r') 



(12) 



In the last line we have given the Fourier decomposition with respect to 
Matsubara frequencies utm — 27Tm/'hf3 omitting the zero mode. Note that 
Gp°'^°{t,t') = —G'p°'^°{t',t). These Green functions possess an interesting 
scaling property: Substituting r = t//3, the Green function becomes indepen- 
dent of 13: 

^p„.o(__,^^_^^2(r-T')-?i[e(T-r')-e(r'-r)]}. (13) 

Introducing expectation values as 



I _ _ fl /■''^ 

2Trh d) VxVp 5{xQ—x)5{pf)—p) ■ ■ -exp \t I ' 



dTip{T)—x{T) 



(14) 

the two-point functions are obtained from functional (jlOp by performing ap- 
propriate functional derivatives with respect to j(r) and v{t), respectively: 

(x(r)i(r'))r^° = 0, (15) 

(i(r)p(r'))r'°=GP»"«(T,r'), (16) 

iPirMrX"^'' = 0. (17) 

The result is that only mixed position-momentum correlations do not vanish. 
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4 Perturbative Expansion for Effective Classical Hamiltonian 

Expanding the restricted partition function ([8]) about a vanishing Hamilto- 



nian, 



X (H(p(T,),i(n)) • • - H(p(t„), i(t,.)))S"'" , (18) 

rewriting this into a cumulant expansion, and utilizing the relation (|5]) between 
restricted partition function and effective classical ffamiltonian, we obtain 

X (iJ(p(ri), x(ri)) . . . Hipir^), a:(r„)))g°;° . (19) 

Using Wick's rule, all correlation functions can be expressed in terms of prod- 
ucts of two-point functions. Since only mixed two-point functions (J12p can 
lead to nonvanishing contributions to the effective classical Hamiltonian, we 
use the rescaled version (fT3|) of the Green function. The scaling transforma- 
tion gives a factor /3 from each of the n integral measures. Thus the expansion 
(|19p is a high-temperature expansion of the effective classical Hamiltonian: 

°° ('—1')"+"'" f^ f^ 

X {H{p{T^), x{t^)) ■ ■ ■ Hip{T„),x{TrMT ' (^0) 

For the following considerations it is useful to assume the Hamilton function 
to be of standard form: H{p{t), x{t)) — p'^{t)/2M + gV{x{T)). Here we have 
introduced the coupling constant g of the potential. Defining the function- 
als a[p] = /p dTp'^{T)/2M and b[x] = J^ dTV{x{T)), the high-temperature 
expansion ((2U)) is expressed as 

n=l ■ k=0 ^ ^ 

In the sequel we point out how this high-temperature expansion is connected 
with an expansion in powers of the coupling constant g of the potential. 
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5 High- Temperature Versus Weak- Coupling Expansion 

In the preceding section it was sliown tliat tlic perturbative expansion about 
a vanisliing Hamiltonian leads to a perturbative series in powers of tire in- 
verse temperature in a natural manner. Now we elaborate the relation to a 
perturbative expansion in powers of the coupling constant g of the potential. 
Changing the order of summation in Eq. (|2ip , one obtains 

(22) 
which is rewritten after explicitly evaluating the (n — 0)- and (fc — 0)- 
contributions: 



2M ^ pf^^^ ^^ n!fc!;i"+'=(2M)» 

X dri • • ■ / dTk dTk+i ■■■ dTk+n 

JQ JQ JO JO 

X (y(a;(ri)) • • • F(a;(T,))p2(r,+i) . • .p2(^,,^„))^°;" (23) 

In this expression, we have inversed the scaling transformation and used that 
j;;^dr {pHr)y;:'' = nppl and /^'^ rfr (y(a:(r)))j;;>- = hl3V{^o). Higher- 
order expectations of functions only depending on a; or p are zero. This is 
due to the vanishing of expectations of functions of x or p which are decom- 
posed into products of two-point functions p3|) and p7)) . All other possible 
contributions are non-connected. 

Note that the expansion ([23)) is equal to the expansion about the free 
particle 

H..(j>o,^o) = ^+9VM + -Y./'-^l d-^-l dr, 

>^{V{x{T,))--.V{x{Tk)))Z,^^. (24) 

where the new cumulants are formed from expectation values 

(• • •)frcc = 27r?i j VxVpSixo - x)S{po -p)--- 



X exp 




d 1 

'i{p{T) -Po)^{x{t) - xo) + —{p{t) -pof 



(25) 
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6 Effective Classical Hamiltonian of Harmonic Oscillator 

In this section, we calculate the effective classical Hamiltonian for the har- 
monic oscillator 

^-(P,^) = |^ + ^M^'^' (26) 

by an exact resummation of the high-temperature expansion (j2ip . For sys- 
tematically expressing the terms of this expansion, it is useful to introduce 
the following Feynman rules: 

ri-^— 72 ^(p(riMr2))r"=Po, (27) 

ri T2 ^(x(Ti)x(T2))r"=^o, (28) 

Ti -— T2 = {x{t{)p{t2))Z' = G''«"«(n,^2) + XmU (29) 

Ti — - T2 =(p(Ti)x(T2))P""'' = -GP°^"(ri,r2)+poa;o, (30) 

r-^ ^(p(r))r"=Po, (31) 

r^ ^(x(r))r"=^o, (32) 

h 
df, (33) 



where the current-like expectations in (|3T|) and (|32|) arise from {p{t))q''^° = 
and {S:{t))q°'^° = 0, respectively. In order to simplify the investigation of the 
expectation values in the high-temperature expansion of the effective classical 
Hamiltonian (1201) . we also define operational subgraphs which are useful to 
systematically construct Feynman diagrams; 



^f^/(^), (34) 



2Mn 



1 /■'* 

— Mw^ / dfx^ir). (35) 

2'i Jo 



Feynman diagrams are built up by attaching the legs of such subgraphs to 
others or by connecting a leg with a suitable current. Note that only com- 
binations of different types of subgraphs lead to non-vanishing contributions, 
since the connection of subgraphs of same type, 

~^.^^^, ^^, (36) 

sets up a new subgraph which contains a propagator (|27|) or ((28|). respectively. 
These propagators, however, are independent of t, such that the r-integrals 
related to the vertices in these subgraphs are trivial. Thus there does not 
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really exist a connection between these vertices and the propagators (^7)1 and 
(p8| can be expressed by the currents (|3T|) and p2|): 



Ti r2 = Ti —* *-~ r2 , (37) 

Ti r2 = ri — * *— r2 . (38) 

Therefore, for n > 1, connected diagrams containing propagators of type (|27p 
or (j28p mwsi break into non-connected parts which we are not interested in. 
Analytically, this is seen by considering for example 

(a;(ri)x(r2))^°^° = {x{T^)iir,)}l°^'> + (x(ri))^«^'' (^(^2))^"" • (39) 

The first term on the right-hand side vanishes due to Eq. (jlSp . while the second 
simply yields Xq, which proves Eq. pop . This means that only Feynman 
diagrams which consist of a mixture of subgraphs (j34p and (j35p contribute to 
the effective classical Hamiltonian. To illustrate this, we discuss the first and 
second order of expansion (|20l) in more detail. 

The Feynman diagrams of the first-order contribution to the effective 
classical Hamiltonian are simply constructed from the subgraphs 

^i!off(Po,a;o) oc — - + ^— 

= 2M,Q ^^'■"■''O '2k'"" ^h"""^~^ 

-ll + iM.^.2, (40) 

where we have used the identities ([57]) and ([55]) in the second expression of 
the second line. Note that the first order term (HUl) obviously reproduces 
the classical Hamiltonian. This is the consequence of the high-temperature 
expansion pip , since only the first-order contribution is nonzero in the limit 
/3 = l/ksT — !■ 0. The second-order contribution reads 

<cff(PO,Xo) ^{--^ +^)(._ +^) 
2 

= — %-{i *-H— ^ -t-4 O- V (41) 

The chain diagram is zero, while the loop diagram has the value — ?i''C(2)/27r^, 
where 



C(-) -f.^ (42) 



n 

ri=l 
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is the C-function. Thus we obtain HJ;^lg{po,xo) = f3h'^uj'^({2)/ATr'^. This 
second-order contribution (|4ip shows the characteristic types of Feynman di- 
agrams in each order n > 1 of the expansion ()20p for the harmonic osciUator: 
chain and loop diagrams. In order to calculate the nth-order contribution, 
we must evaluate these diagrams more general. By constructing Feynman 
diagrams from the product of n sums of subgraphs, 



)( 



)•••( 



+ 



n times 



it turns out that only following chain and loop diagrams contribute: 



(43) 



\ 



(44) 



The evaluation of the chain diagrams is easily done and yields zero. An explicit 
calculation in Fourier space shows that there occur Kronecker-(5's 6m o- Since 
the Matsubara sum of the Green function (HH)) does not contain the zero 
mode m — 0, all chain diagrams are zero. 

Determining the values of loop diagrams is more involved. It is obvious 
that loop diagrams can only be constructed in even order {n — 2,4, 6, . . .), 
since for a loop diagram with mixed propagators (l29l) or (l30l) pairs of different 
subgraphs (p4)) and (|35t are necessary. Thus we have found the result that odd 
orders of expansion (|2T|) vanish, and only loop diagrams for n G {2,4,6,...} 
must be calculated. Evaluating loop diagrams of nth order in Fourier space 
is straightforward and entails 



2(-l) 



2fc 



^) «»>• 



(45) 



where k = n/2. The high-temperature expansion for the effective Hamiltonian 
of the harmonic oscillator can thus be written as 



pI 



Hu:,cff{po,Xo) ^ f^ + l^Mu?xl 



fc=l 



p' 



(-1) 



fc+1 



2k 



g) C(2fc). (46) 



Substituting the (^-function by its definition (|42p and exchanging the summa- 
tions, the last term in Eq. (|46l) can be expressed as a logarithm 

k V 27r 



Y^ o2fc-l 



fc=l 



/3 



C(2fc) = -In 



n 



2«2, ,2 



1 



h'(3^u; 



Air'^n? 



(47) 
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Applying the relation 

-sinh.= n l + -i-I , (48) 



n=l ^ '' 

we find the more familiar form of the effective classical Hamiltonian for a 
harmonic oscillator 

«„.»(po,.o)^|| + iM.^.2-i,„^-5^, (49) 

Performing the xq- and po-integrations in Eq. ([T]), we obtain the well-known 
form of the partition function of the harmonic oscillator Z^^ — 1/ sinh huj/3. 

7 Summary 

We have used a perturbative definition of the path integral in phase space 
representation which reproduces the effective classical Hamiltonian for the 
harmonic oscillator. Our procedure is an alternative way to evaluate path 
integrals: The unperturbed system is trivial and the calculation of appro- 
priate Feynman diagrams is simple. Furthermore it turns out that the per- 
turbative expansion for the effective classical Hamiltonian is identical to the 
high-temperature expansion. 
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